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Abstract 

Numerical solutions of the compressible 
Navier-Stokes equations by an Alternating Direc- 
tion Implicit scheme, applied to two experimen- 
tal investigations at NASA Langley Research 
Center are presented. The first is cooling by 
injection of a gas jet through the nose of an 
ogive-cone, and the second is the aerothermal 
environment in the gap formed by the wing and 
elevon section of a test model of the Space 
Shuttle. The simulations demonstrate that 
accurate pressure calculations are easily 
obtained on a coarse grid, with convergence 
being obtained after the residual reduces by 
four orders of magnitude. However, accurate 
heating rates require a fine grid solution, with 
convergence requiring at least a reduction of 
six orders of magnitude in the residual. The 
effect of artificial dissipation on numerical 
results is also assessed. 

I. Introduction 

In recent years an increased interest in 
hypersonic vehicle concepts has brought a number 
of aerodynamic problems for hypersonic flow into 
prominence. These problems are characterized by 
high-temperature high-velocity flows with strong 
viscous interactions such as shock waves, blunt 
leading edges, etc. A dominant feature of these 
flows is the aerodynamically induced high heat 
transfer rates to the vehicle surfaces. Accu- 
rate heat transfer rates are difficult to obtain 
numerically, and difficult to measure experimen- 
tally. However, accurate heating rates are 
important when designing a given structure and 
assessing its reliability. 

Numerical and experimental models are 
important tools for understanding the aerother- 
mal environment in high-speed flows. However, 
before extrapolating the results of such models 
to actual flight conditions careful model vali- 
dation along with an understanding of the 
effects model assumptions have on the simulated 
physics are important first steps. 

This paper will present results from numer- 
ical models of two experiments investigated at 
NASA Langley Research Center. 1 The studies 
were undertaken to complement experimental 
results by providing detailed information not 
obtainable experimentally. In addition, these 
studies assess the effects artificial dissipa- 
tion, convergence criteria, and degree of reso- 
lution have on pressure and heating rate predic- 
tions, and hence provide guidelines for improved 
techniques. The solution technique is a 
second-order accurate Beam and Warming 
Alternating Direction Implicit (ADI) 
algorithm. 2 The algorithm is augmented to 
include an explicit blend of second and fourth 
differences for the dissipative terms.- 1 
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The first problem is Mach 6.7 flow over a 
gas jet cooled ogive-cone (gas jet nose tip, 

GJNT ; see Fig. la). Protecting a nose tip from 
severe thermal and pressure environments is an 
important concern for advanced entry vehicles. 
Very little experimental or numerical results 
exist for high enthalpy hypersonic flow condi- 
tions for GJNT on large cones. Therefore, a 
test program was conductd in the Langley 8-Foot 
High Temperature Tunnel (8‘ HTT) to study mass 
addition cooling by GJNT.4»5 Numerical shock 
stand-off distances, pressures and cold wall 
heating rates will be compared to experimental 
results. 

The second problem studied in only a preli- 
minary fashion, is the determination of the flu- 
id/thermal environment in the cove between the 
Space Shuttle wing and elevon surfaces (Fig. lb) 
which has been a concern throughout the design 
and initial operational phases of the vehicle. 

If the seals at the elevon hinge line leak, hot 
boundary layer gas flowing into the cove will 
endanger the thermally unprotected inner wing 
and elevon structures. Previous attempt using 
finite difference techniques 6 were unable to 
obtain quantitative results, possibly due to 
difficulties encountered with modelling the com- 
plex geometry. Finite element techniques have 
handled the geometric complexity and solutions 
have been obtained for both incompressible flow 
in the cove only 7 and for inviscid flow for 
the external and internal flow fields.* 1 How- 
ever, since details of the flow at the mouth and 
internal to the cove will be dominated by vis- 
cous effects and the shear layer emanating from 
the rear of the wing panel a full Navier-Stokes 
solution is required. The numerical analysis to 
be given models the flow with the viscous-com- 
pressible Navier-Stokes equations. Comparison 
of results with data from a full scale model for 
the wing-elevon cove tested in the 8' HTT^ is 
also given. 

II. Numerical Model 

The basic equations governing the flow of a 
viscous compressible fluid are the Navier-Stokes 
equations. ^ These equations, expressed in 
terms of a two-dimensional or axi symmetric curv- 
ilinear coordinate system, are written in strong 
conservative form for a perfect gas as follows: 

U t + h + F n + jH = m (R C + S n + J ’ T)/Re (1) 

with 

U = U/J 

E = Ut H + ^x E - + 

F - (nt U + nx E + OyF)/J 
H = (F +J)/(yJl 7 
R - U x 3 + Sy §)/J 
S = (fiy R + Tiy S)/J 
T = (S + T)/(yJ) 

where U, E, F, H, R, S, and T are the following 
four-component vectors: 
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and t y y = (x + 2 y) u x + xv v + jxv/y 
T Xy =u(u y +v x ) y 

Tyy = (X + 2 \l) Vy + XU X + jxv/y 

T 0 0 = (x + 2 u) v/y + X (u x + Vy) 


In these equations, £ and n are the computation- 
al coordinates, t is the time and J is the Jaco- 
bian of the coordinate transformation. The com- 
putational coordinates are related to Cartesian 
coordinates through the coordinate transforma- 
tion 


5 = 5 (t,x,y) (2) 

n = n (t,x,y) 

(For generality the coordinates are written as a 
function of time though the cases which follow 
are steady-state.) In Eq. (1) j=0 for two- di- 
mensional flow, j=l for axi symmetric flow (y be- 
comes the cylindrical radius), m=0 for inviscid 
flow, and m=l for viscous flow. 

In the conservative variables of Eq. (1), 
density, velocities, and temperature are nondi- 
mensional ized with respect to freestream 
values. The pressure, density, and velocity 
components are related to the total energy per 
unit volume, e, by the following equation for an 
ideal gas: 

e=p [e+(u 2 + v 2 )/2] = p/(y-l) + p(u 2 + v 2 )/2 (3) 

Re represents the Reynolds number. The Prandtl 
number, Pr, the coefficient of thermal conduc- 
tivity, k, and the ratio of specific heats, y , 
are held constant. The coefficient of thermal 
conductivity and the viscosity coefficient are 
nondimensional ized with respect to their free- 
stream values. The viscosity coefficients, y 
and X, are computed from Sutherland's equation 
and Stokes' hypothesis, respectively. ^ Boun- 
dary conditions will be discussed when a given 
problem is addressed. 

The solution technique is a second order 
accurate Beam and Warming ADI algorithm. The 
algorithm is augmented to include an explicit 
blend of second and fourth differences for the 
dissipative terms regulated by the pressure gra- 
dient in the local flow field. Their addition 
controls the odd-even uncoupling of grid points 
typical of central difference schemes, and pro- 
vides strong dissipation in the region of shock 
waves and lesser dissipation in the remainder of 
the field. 


All calculations were performed on the 
Langley Cyber 175 and took ,00247sec/pt ./itera 
tion of CPU time. 

III. Gas Jet Nose-Tip 


Nonbl owing Cases 

Before modelling the GJNT, a series of cas- 
es for cones without gas injection were studied 
to assess the effects that convergence, resolu- 
tion, and artificial dissipation have on numeri- 
cal results without the added complications of 
the gas jet. Freestream conditions are M m = 

6.7, Re„=1.4xl0 6 ft.' 1 , and T„ = 407°R. 

Boundary conditions are no-slip at solid sur- 
faces, extrapolation at outflow and a wall temp- 
erature of 540°R. 

Numerical results are compared with experi- 
mental data obtained in the Langley Research 
Center 8-Foot High Temperature Tunnel. 12 The 
facility is a hypersonic blowdown wind tunnel 
that operates at a nominal Mach number of 7, in 
a test medium of combustion products of methane 
and air. Numerical results will be compared 
with two nose configurations on a 12.5° half- 
angle cone: a 3 inch radius tip (R3) and a one 

inch radius tip on an ogive frustrum (Rl). Un- 
less otherwise stated the following results and 
comparisons are for Rl, the GJNT configuration 
without the jet. 

Fig. 2 is a plot of shock stand-off dis- 
tance comparing numerical results with experi- 
mental data. 4 As indicated in Fig. 2, stand- 
off is measured as the distance d^, to the 
body divided into the distance d 2 to the bow 
shock. Distances are measured from the symmetry 
line of the body. Numerical standoff measure- 
ments are about 7% higher than experimental mea- 
surements. The impact of grid resolution on 
capturing this shock is illustrated in Figs. 3a 
and 3b by Mach numbers contours. The fine grid 
solution has 35 points along the body and 41 
normal to the body; and the coarse grid (17 x 
21) is obtained by deleting alternating points 
from the fine grid. The shock tends to steepen 
up towards the low Mach number side which exper- 
ience has shown to be true of central difference 
schemes. The shock width spans three grid cells 
in both grids. It is interesting to note a re- 
sult correlating numerical results for different 
gases and hemispherical bodies. 13 The corre- 
lation is a simple linear law relating shock 
standoff distance to the density ratio across a 
normal shock. Fig. 3c is a magnification of 
Fig. 3b illustrating the shock location from 
reference 13 and the experimental shock location 
at the stagnation point. The correlation pre- 
dicts the shock to lie on the low Mach number 
side of the present numerical model. It should 
be noted that the numerical model of reference 
13 was inviscid and adiabatic, which indicates 
that the shock stand-off is not a strong func- 
tion of these factors. 

Fig. 4 is a plot of calculated surface 
pressures from the fine grid after a four order 
of magnitude reduction in the residual, compared 
with experimental data. The abscissa is surface 
distance measured in units of nose radii. Pres- 
sures are nondimensional ized by stagnation point 
values. The ratio of stagnation pressure to 
freestream static pressure Is 60.67 for the nu- 
merical model and 61.90 for the experiment. The 
experimental and analytical data overlap between 
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three and ten inches and both are seen to fall 
on the same curve. 

The effect of convergence criteria on 
pressure is illustrated in Fig. 5 which is a 
plot of pressure profiles from runs using coarse 
grids and terminated at three different average 
global residuals. Stagnation pressures divided 
by the freestream pressure for residual levels 
of 10* 4 , lO- 6 , and 10* 7 are 57.89, 56.78 
and 56.77, respectively. A four order of 
magnitude reduction in the residual (2000 
iterations) differs very little from six and 
seven orders of magnitude reductions. 

Converged solutions of pressure are used to 
illustrate effects of grid resolution in Fig. 

6. The ratio of stagnation pressure over free- 
stream pressure is 56.78 for the coarse grid 
solution and 60.67 for the fine grid case. Dif- 
ferences in the two solutions for pressure are 
small, the discrepancy between 1 1/2 and 2 nose 
radii reflect the inadequacy of the coarser grid 
in the tangential direction to resolve the ex- 
panding flow around the hemispherical nose. 

Fig. 7 is a plot of numerically predicted 
pressures (residual = 10” 6 ) from a coarse grid 
solution versus surface distance for different 
levels of artificial dissipation. The effect of 
doubling the dissipation has a noticeable effect 
on the pressure ratio profiles: higher dissipa- 

tion resulting in higher pressure profiles. 

This effect illustrates the impact artificial 
viscosity has on the inviscid portion of the 
flow field where second order dissipation terms 
are needed to capture the bow shock. The magni- 
tude of the dissipative terms should be kept at 
the minimum required to control oscillations and 
maintain stability of the scheme. However, de- 
termining the amount required for a particular 
application is typically a matter of exper- 
ience^ and dependent on the configuration, 
solution algorithm and dissipation model. The 
value of the stagnation pressure divided by the 
freestream pressure is 56.78 and 53.09 for the 
low and high dissipation cases, respectively. 

A comparison of calculated heating rates 
and experimental data, each normalized by its 
own stagnation point value for R3, is shown in 
Fig. 8. Since the governing equations are in 
nondimensional form, the ratio of heat transfers 
is the ratio between nondimensional temperature 
gradients. Experimental heating rates are non- 
dimensional ized by stagnation point values cal- 
culated from reference 15. The wall temperature 
for the Navier-Stokes calculation is fixed at 
540°R. This is the wall temperature at which 
the experimental data was extrapolated as a 
standard reduction procedure. Experimental data 
begin approximately three nose radii from the 
tip. Numerical results are seen to overlap the 
experimental data in this region. Fig. 9 is a 
similar plot for Rl. 

Convergence criteria are found to be more 
strict for heating rates than for pressures. At 
least a six order of magnitude reduction (5000 
iterations) in the average global residual is 
needed to converge heating rates to a steady- 
state value. Fig. 10 shows heating rates con- 
verged four, six, and seven orders of magni- 
tude. Seven orders of magnitude reduction shows 
almost no change from the six order of magnitude 
reduction. The normalization factor for all 
three runs are within 1.4% of each other. It 
should be noted that the Beam and Warming 
algorithm performs faster for problems with 


weaker shocks, though typically an order of 
magnitude slower than present day upwind 
schemes. 16 Upwind schemes, however, have not 
yet been tested on complicated flow fields. In 
addition recent work by Pulliam, 14 indicates 
an order of magnitude increase in convergence 
rate for the Beam and Warming algorithm with the 
addition of implicit fourth order dissipation. 

The effect of grid resolution on a fine 
(35x41) and coarse (17x21) grid for converged 
heating rates is illustrated in Fig. 11. The 
normalization factor is 10% lower for the coarse 
grid solution. The disagreement is greatest 
between one and two nose radii from the tip. 

This reflects the inadequacy of the coarser grid 
in the tangential direction to resolve the flow 
expanding around the hemispherical nose. 

Fig. 12 illustrates the effect doubling the 
dissipation has on heating rates for coarse grid 
solutions which are fully converged (residual = 
10“'). The profiles are not identical; the 
higher dissipation resulted in lower heating 
rates. The lower dissipation case has a normal- 
ization factor which is 7 % higher than the high 
dissipation case. 

Gas Jet Nose-Tip 

To reduce the excessive heating loads on 
bodies in hypersonic flow a study of the cooling 
effectiveness produced by mass injection was 
conducted in the 8' HTT.^» 5 A parallel study 
is presented using the Navier-Stokes code. 
Freestream conditions are H* « 6.7, R^ = 1.4 
x 10 6 ft.* 1 , y=1.382, Prandtl number =.743. 

The gas jet has a temperature of 495°R and is 
sonic at the outlet. The total pressure of the 
coolant divided by the pitot pressure of the 
flow is 2.46 producing an underexpanded jet of 
coolant gas. These conditions are the same as 
tunnel test conditions. However the experiment 
is in a test medium of combustion products of 
methane and air and the gas jet is nitrogen, 
while the numerical model assumes perfect gas 
relations and ignores molecular diffusion. 
Numerical results are obtained on a 35 x 41 grid 
(35 points along the body and 41 normal to the 
body) and residuals are reduced by 6 orders of 
magnitude. 

In 1965 GJNT experiments conducted by 
Finley 17 on a blunter body led him to postu- 
late flow features based on Schlieren, shadow- 
graph and experimental data. A schematic taken 
from ref. 17 is reproduced in Fig. 13. The 
schematic depicts a low velocity recirculation 
zone surrounded by a shear or mixing layer which 
attaches to the body beyond the recircul ation 
zone. A comparison between GJNT experiments 
conducted in the 8' HTT and the numerical model 
reveals similar features. 

Fig. 14 is a qualitative comparison of 
shock stand-off distance by interfacing a 
shadowgraph from the GJNT experiment 2 with 
Mach number contours from the numerical solu- 
tion. Stand-off distance compares well with the 
low Mach number side of the simulated shock. 
(Recall from Fig. 3 that a finer grid will cause 
the shock to steepen up towards the low Mach 
number side.) Three cells span the shock width 
which is expected with first order accuracy in 
this region. A quantitative plot of shock 
stand-off distance versus angular distance about 
the symmetry line at the nozzle exit plane is 
given in Fig. 15. 
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Fig. 16 is a plot of temperature contours 
(T w = 540°R) which helps illustrate some 
important flow features. The distance indicated 
as s=4.0 inches (axial) encompasses a recircula- 
tion region, as confirmed by numerical values of 
velocities, where the cold gas is most effective 
in cooling the body. Downstream of this recir- 
culation zone, at approximately four inches from 
the nose tip, flow reattaches and the tempera- 
ture contours move in toward the body though a 
cold layer of gas still persists down the length 
of the body. Heating rates remain negative out- 
side the recirculation region, but increase 
after the recirculation zone. Some of these 
same features are revealed by image enhancements 
of experimental photographs. Fig. 17 is an 
image enhancement of a Schlieren taken of a GJNT 
experiment. 5 The contrasts help delineate the 
recirculation region and an outer mixing layer 
of hot incoming gas mixing with the jet and 
impinging the body downstream of this recircula- 
tion zone. Fig. 18 is a composite of two 
enhancements of the shadowgraph in Fig. 14. It 
shows further detail of the recirculating gas 
jet and depicts a turbulent mixing layer outside 
the recirculation region. Experimental heat 
transfer data also indicates turbulent mixing to 
be occurring. These observations are not sur- 
prising since the Reynolds number of the gas jet 
based on the nozzle diameter is approximately 2 
x 10°, which is three orders of magnitude 
higher than the Reynolds number for fully- 
developed turbulent pipe flow. Fig. 19 is a 
comparison of numerical and experimental heating 
rates. These heating rates are nondimensiona- 
lized by stagnation point values of a nonblowing 
experimental case. The heating in the numerical 
model is due to lateral conduction as expected 
in a laminar flow. The experimental data show 
high positive heating rates downstream of the 
recirculation zone indicating turbulent mixing 
of the hot incoming gas with the turbulent gas 
jet. Note that experimental heating rates 
radically increase at approximately four inches 
from the nose tip, which is the length of the 
recirculation region indicated in the tempera- 
ture contour plot of Fig. 16. However, the 
numerical heating rates remain negative down- 
stream of the recirculation zone. The numerical 
results for the adiabatic wall case show wall 
temperatures below 540°R which confirm the nega- 
tive heating rates in Fig. 19. Fig. 20 is a 
plot of adiabatic wall temperatures for a cold 
as jet (495°R) and a slightly warmer gas 
540°R) . The adiabatic temperatures are well 
below the fixed cold wall temperature, the 
warmer gas having slightly higher adiabatic 
temperatures. As expected a warmer gas results 
In higher heating rates as shown in Fig. 21, 
which plots heating rates for gas jets of 495°R 
and 800°R . Figs. 20 and 21 reveal that differ- 
ences in gas jet temperatures are not sufficient 
to account for the discrepancy between numerical 
and experimental heat transfer since the 
Increase in heating rates is small for a 395° 
rise in gas jet temperatures . Thus, the numer- 
ical comparison with experiment suggests that 
maintaining the gas jet as laminar as possible 
will increase its cooling effectiveness. 

Comparison of experimental ^ and numerical 
pressures normalized by stagnation point values 
of a nonblowing case is given in Fig. 22. The 


numerical pressures indicate a dip in the 
recirculating region. Experimental data begin 
outside of this region at approximately 12 
inches. The comparison indicates good agree- 
ment. A calculation of the ratio between the 
pressure behind the normal shock and P M for a 
nonblowing case is 57.702. The highest pressure 
ratio in the GJNT numerical model, experienced 
at the lip of the ogive, is around 31. Thus the 
gas jet substantially reduces the pressure load 
on the nose tip. 

IV. Preliminary Study of a Wing-Elevon Cove 

The last problem to be addressed is a study 
of the fluid/thermal environment in the cove 
between the wing and elevon on the Space Shuttle 
(Fig. lb). The gap between the Space Shuttle 
wing and elevon surface is closed by seals at 
the elevon hinge line to prevent the leakage of 
hot boundary layer gases. If these hot gases 
flow into the cove, damage to the thermally 
unprotected inner wing and elevon structures 
could occur. To provide insight into the prob- 
lem a full scale model of the wing-elevon cove 
was tested in the 8‘ HTT.^ A 2D Navier-Stokes 
solution of the aerothermal environment of the 
wing elevon cove is given. A major difficulty 
in obtaining a reliable numerical solution has 
been producing a computational grid which models 
the cove geometry adequately. A fully interac- 
tive 2D algebraic grid generator 18 currently 
under development by Erlebacher at Langley, is 
used to generate the grid required for the 
present work. Fig. 23a show the grid produced 
by this method, for a cove geometry in which the 
seal "leak" area is equal to the entrance area 
and the elevon deflection is 25° with respect to 
the wing. The grid is composed of 55 points on 
the coordinate family which spans the cove gap, 
and 81 points on the opposite family. The 
entrance "neck" region is quite difficult to 
grid with adequate resolution and near- 
orthogonal ity ; the enlargement of that region 
shown in Fig. 23b shows good performance by the 
above procedure. 

In the experimental investigation of the 
configuration, the elevon/gap apparatus was 
mounted in a flat test bed inserted into the 
tunnel; a flat run of approximately 48 inches 
preceded the cove entrance. The entire test bed 
was inclined 5° (compression) with respect to 
freestream. The tunnel freestream conditions 
were = 6.7, R*, = 1.0 x 10 8 ft -1 and 
^ = 399°R; the Reynolds number based on cove 
entrance gap was 115,000. 

Upstream inflow conditions were difficult 
to set in the numerical model, since the 
inclination of the experimental apparatus 
resulted in a weak shock, with a slight entropy 
layer due to the small radius (3/8") leading 
edge. Since no detailed surveys of the flow- 
field upstream of the cove were available, 
numerical inflow conditions were approximated by 
setting tunnel freestream conditions inclined at 
5° to the wing surface, blending into a flat- 
plate boundary layer profile near the wing sur- 
face. The thickness of the boundary layer was 
about 20% of the cove entrance gap, consistent 
with the flat plate run preceding the cove 
entrance in the experimental apparatus. The 
effect of these approximate inflow conditions 
requires further assessment. No-slip, flow 
tangency and zero normal -pressure gradient were 
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imposed on the wing and elevon surfaces, along 
with either adiabatic (zero normal -temperature 
gradient) or constant wall temperature (540°R) 
conditions. Extrapolation was used at both 
freest ream outflow and cove "leak*' regions. 

Fig. 24 is a Mach number contour plot of the 
flow in the wing-elevon cove. Note the weak com- 
pression waves upstream of the cove entrance due 
to the 5° angle of the input conditions as men- 
tioned previously. Since in the experiment the 
leading edge is 48 inches from the cove entrance 
the compression waves did not have a comparable 
impact on the flow in the near vicinity ofthe 
cove entrance. The boundary layer separates from 
the wing trailing edge forming a shear layer 
across the cove entrance which impinges on the 
deflected elevon downstream where the elevon 
shock forms. Separation occurs just upstream of 
the cluster of flow direction arrows at the wing 
trail ing-edge in Fig. 25 indicating a reversal 
of flow direction. On the cove side of the 
shear layer a large primary recirculation region 
forms and flow is directed into the cove. A 
small secondary recirculating area appears on 
the elevon surface inside of the cove. In Fig. 
25 this secondary vortex is seen as the flow 
initially heads up the cove and then abruptly 
reverses direction in a confined region above 
the recirculation. The primary recirculation is 
quite strong; local Mach numbers exceed one in 
some areas. However, little flow is induced 
farther up in the cove, where Mach numbers are 
less than .3. 

Comparisons are made with experimental data 
in Fig. 26 for the pressure coefficient C p 
versus surface distance. As indicated in rig. 
23a measurements begin at 0 in the cove and pro- 
ceed down the elevon surface. Numerical values 
of Cp are higher than experimental data. The 
prediction of the slope for the steep C p rise 
is higher than the experiment. This rise occurs 
where the shear layer merges with the shock. 

The occurrence of transition in this separated 
shear layer is likely to influence the value of 
this gradient. Experimental results for this 
test case 9 indicate heat transfer data on the 
elevon to be within a factor of two of fully 
turbulent flow as calculated by a reference en- 
thalpy method, thus suggesting the likelihood 
that this separated shear layer is indeed in 
transition. It should be noted that a similar 
effect might occur if the experimental boundary 
layer was substantially thicker than the one 
assumed in this analysis, (the lower velocities 
causing lower pressure recovery on the elevon 
and a longer distance before the inviscid level 
is obtained.) 

The value of the C p (i.e. pressure) pla- 
teau on the elevon is within a few percent of 
the experimental value verifying proper post 
shock inviscid conditions in the numerical model 
The dip in Cp at approximately 7.5 Inches is 
at the attachment point. There are not enough 
experimental data points to plot a profile in 
this section. 

VII, Conclusions 

Two aerothermal flow problems of relevance 
in hypersonic flow research have been modelled 
by compressible Navier-Stokes calculations using 
an ADI scheme and compared with experiment. In 
this process, the impact of resolution, conver- 
gence and dissipation on predicted pressures and 


heating rates was assessed. Heating rates were 
found to require at least a six order of magni- 
tude reduction in residual for full convergence 
whereas pressures were converged after a reduc- 
tion of four orders of magnitude. In addition, 
accurate heating rate predictions required a 
finer grid than accurate pressure predictions. 

In both cases, discrepancies in the fine and 
coarse grid solution occurred at the stagnation 
point and where the flow expanded rapidly, re- 
flecting the inadequacy of the coarse grid in 
the tangential direction to resolve the expan- 
sion. Doubling both the second and fourth order 
dissipation produced noticeable changes in both 
pressures and heating rates. The higher dissi- 
pation reduced predicted heat transfer rates and 
increased pressures. Such effects illustrate 
the impact artificial dissipation has not only 
on the boundary layer but in the inviscid por- 
tion of the flow field where the second order 
dissipation terms are needed to capture shocks. 

The ability of the Navier-Stokes/ADI scheme 
to model complex hypersonic flows was illustrat- 
ed by comparison with experiment. A numerical 
model of the GJNT at zero angle of attack com- 
pared favorably with experimental pressures, 
shock stand-off and qualitative features. How- 
ever, experimental heating rates beyond the re- 
circulation zone were much higher than the nu- 
merical heating rates. Differences in the gas 
jet temperature was shown to be insufficient in 
the numerical model to account for the large 
discrepancy in the comparison. However, thenu- 
merical model is laminar and the heating rates 
from the experiment may be dominated by turbu- 
lent mixing of the incoming hot gas with the 
turbulent gas jet. The considerably lower heat- 
ing rates in the numerical model suggests that 
maintaining the jet as laminar as possible might 
greatly enhance its cooling effectiveness. 

Preliminary results from numerically model- 
ling a wing-elevon cove indicate good agreement 
with experimental values of C p and qualitative 
flow features. An over prediction by the numer- 
ical model in the magnitude of a steep C p gra- 
dient occurring in the vicinity of the separ- 
ated shear layer may be due to transition ef- 
fects in the experimental model, since experi- 
mental results indicate heat transfer data on 
the elevon to be within a factor of two of fully 
turbulent flow. In addition numerical inflow 
conditions were approximations to the conditions 
present in the experimental model; the effect of 
proper inflow conditions requires further as- 
sessment. However, the value of the C p pla- 
teau after its steep rise, predicted by the nu- 
merical model is within a few percent of the ex- 
perimental value verifying the proper post-shock 
inviscid conditions in the numerical model. 
Agreement with experimental data in both cases 
would be greatly enhanced by including turbu- 
lence and real gas effects in the numerical 
model . 

For complicated flow fields Beam and Warm- 
ing is a proven algorithm in performance and re- 
liability. Present day upwind schemes have not 
yet been tested for complex flow fields, though 
for simpler problems upwind schemes do exhibit 
an order of magnitude faster convergence rate 
than Beam and Warming. However, recent work by 
Pulliam indicates an order of magnitude increase 
in speed after the addition of implicit fourth 
order dissipation to the Beam and Warming algor- 
ithm, which would make it competitive with 
upwind methods. 
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Fig. 1. Tests in the 8-Foot High Temperature 

tunnel at NASA Langley Research Center 

a) gas jet nose-tip ogive-cone 

b) wing-elevon cove 
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Fig. 2. Shock Stand-Off Comparison on 

Nonblowing Ogive-Cone (Rl): Numerical 

and Experimental Results 
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Fig. 4. Calculated and Experimental Pressures 
for Nonblowing Case (Rl) 
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3. Coarse and Fine Grid Solution for Rl: 
Comparison of Shock Stand-Off Distance 

a) coarse grid solution 

b) fine grid solution 

c) magnification of fine grid solution 
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Fig. 5. Effect of Convergence on Pressures (Rl) 
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Fig. 6. Coarse and Fine Grid Solution of 
Pressures (Rl) 
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Fig. 7. Effects of Dissipation on Pressure (Rl) 


Fig. 10. Effect of Convergence on Heating Rates 
(Rl) 
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Fig. 8. Calculated and Experimental Heating 
Rates for Nonblowing Cone (R3) 
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Fig. 9. Calculated and Experimental Heating 
Rates for Nonblowing Ogive-Cone (Rl) 


Fig. 11. Coarse and Fine Grid Solution of 
Heating Rates (Rl) 
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Fig. 12. Effect of Dissipation on Heating Rates 

(Rl) 
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Fig. 22. Calculated and Experimental Pressures 
for GJNT 



Fig. 24. Mach Number Contours of the Wing-El evon 
Cove 
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Fig. 23. Grid System for Wing-Elevon Cove 
Problem 

a) entire computational domain 

b) expanded portion of the "neck" 
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Fig. 25. Flow Direction Plot in Expanded Portion 
of the "Neck" Region 
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